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^ 1 Abstract 

o , 

£S) , It is shown that a new light gauge boson X which might be produced in the decays 

of pseudoscalar mesons 7r (77) — > 7 + X could be effectively searched for in neutrino 
experiments via the Primakoff effect, in the process of X + Z — ► vr°(r/) + Z conversion 
I/"") , in the external Coulomb field of a nucleus. An estimate of the X — ► 7r° conversion rate 

for the NOMAD neutrino detector at CERN is given. 

m ■ 

O . 

00 • 1 Introduction 

ON 

New neutral gauge bosons X are predicted by many models addressing the physics beyond 
the Standard Model (SM) such as GUTs 0, supersymmetric @, superstring models |§ and 
models including a new long-range interaction, i.e. the fifth force, Q]. The predictions for 
■ the mass of the X boson are not very firm and it could be light enough (Mx <C Mz) for 
searches at low energies. 

The detailed study of the possible manifestations of light gauge boson was performed in 
[H], It was shown that if the mass Mx is of the order of the pion mass, an effective search 
could be performed for this new vector boson in the radiative decays of neutral pseudoscalar 
mesons P — > 7 + X, where P = 7r°, rj, or 77'. 

From the analysis of the data from earlier experiments, constraints on branching ratios 
for the decay of P — > 7 + X range from 10 -7 to 10 -3 depending on whether X interacts with 
both quarks and leptons or only with quarks. In the first case X is a short lived particle 
decaying mainly to e + e" or vv pairs, while in the second case X should be a relatively long 
lived particle (tx > 10~ 6 sec) 0. 

Direct searches for a signal from 7r° — >• 7 + X decay have been performed in a few 
experiments. The best experimental limit on the decay 7r° — >• 7 + X was obtained recently by 
the Crystal Barrel Collaboration at CERN using pp annihilations as a source of pseudoscalar 
mesons, see @ and references there. The branching ratio limit of (3 to 0.6) x 10 -4 (90% C.L.) 
has been obtained for < Mx < 125 MeV. Limit of the same order has been obtained for 
BR(i] — > 7 + X) for the higher Mx mass region. The results are valid for the case where X 
is a long lived particle or the X boson decays preferably into uT7 pairs. 

In this paper we show that a new light relatively long lived gauge boson which might be 
produced in decays of 7r°(j]) — > 7+X could be effectively searched for in neutrino experiments 
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via the Primakoff effect, in the process of XZ — > 7r°(r])Z conversion in the external Coulomb 
field of a nucleus. 

The paper is organised as follows. In section 2 we discuss the X-boson phenomenology. 
Section 3 describes the suggested method and presents the results of calculations of the 
cross section for X — ► 7r° conversion. An estimate of the total cross section for the X- 
boson interaction with the matter is given in section 4. In section 5 we consider NOMAD 
neutrino detector at CERN as example to estimate expected X — > 7r° conversion rate in this 
experiment. Section 6 contains concluding remarks and discussion. 



2 X-boson phenomenology 

As it has been mentioned in the introduction many modern models predict an enlargement 
of the standard SU{3)®SU (2) ® U(l) gauge group by an extra £7(1) factor. At present there 
are no firm theoretical prediction on the mass of new gauge boson so we can't exclude the 
possibility that new gauge boson is rather light with a mass Mx < 0(m n ). The possibility 
of the existence of a vector boson, which interacts with the baryon current, was suggested 
by Lee and Yang ||. If this boson is massless, then the results of the experiments on 
the check of the equivalence principle imply a strong constraint on its coupling constant 
oliy < 1CT 47 . However this constraint is no longer valid for massive gauge boson with 
microscopic Compton wavelength. Astrophysical bound on the X-boson coupling is [|1(| 
«! < O(10~ 10 ) for M x < 0(100) KeV. 

The interaction Lagrangian of the X-boson with quarks and leptons has the form [H] 



J X,ql 



where (gv q , 9A q , gvu 9ai) ~ 0(1). By C-parity arguments, one can find that a coupling of 
X-boson to the axial quark current gives a negligible small contribution to the tt° — > 7 + X 
decay width. The tt° — > 7X decay width is determined by the formula || 

r(7r _ 7 + X) = -^—^(2g Vu + gvd ) — (1 - _) , (2) 
x = 9 i,U = 93MeV andM, 

we have 



where ax = f^, fn = 93MeV and Mx is the mass of the X-boson. For the Br(n° — > 7 + X) 



Br (n° ^ 1 + X) = 2^(2g Vu + g Vd )\l - ^ff (3) 

For the case when X-boson interacts with both quarks and leptons combined bounds 
from anomalous magnetic moment of leptons, elastic v e e — > v e e scattering and beam dump 
experiments lead to 

a x <O(10- 9 ), (4) 

or to 

Br(n° ^ 7 + X) < O(10" 6 ) (5) 

The X-boson lifetime for M x > IMeV and a x = 0(1(T 9 ) is less than r(X) < O(10" 12 ) 
sec. So, X-boson is short lived particles and such scenario is not interesting for us. 

In this letter we consider the case where the X is a "leptophobic" boson which interacts 
only with quarks. Consider the life time of the X-boson for such scenario. Experimental 
limits for the Br(7i°(r]) — > 7 + X) lead to the bound on coupling of X to quarks of 

a x < KT 6 (6) 



For Mx < 2m e the X-boson decay width into 3 photons via quark loop is proportional to 



T(X - 37) ~ a 3 a x (^^) s M x , (7) 
m g,eff 

where ?7i ? , e // ~ 300MeV is an effective quark mass. Using Eq.(6) we find that rx(X — ► 
37) > O(10) sec. The decay width of the X-boson into vv via quark loop is proportional to 

T(X -> j/P) ~ a x G^Mi(^^) 4 (8) 



For M x < 100 MeV we find that r x (X -> z/p) > O(10" 4 ) sec. For M x > 2m e the X-boson 



decays mainly into e + e through the quark loop and its decay width for M\ <C rn 2 „ e ff ls 



T(X -> eV) ~ T ^-( T ^) 2 [3(|^ U - - k) 2 ]VM x (i^) 4 (9) 

1/7T 157T 3 3 3 m q,eff 

For M x < 100 MeV and [3(|c/y u - \g Vd - \g V s) 2 ] 2 = 0(1) we find that r x {X -> e+e") > 
O(10 -6 ) sec. So, we conclude that for the case when the X-boson interacts only with quarks 
it is long lived particle at least for Mx < 100 MeV and for X boson energy Ex > 10 GeV it 
travels the distance more than 10 km before its decay which is much larger than the typical 
decay length (~1 km) in the current neutrino experiments. 



3 Method of search and cross section for X —> tt° con- 
version 

Let us consider for simplicity the decay tt° — ► 7 + X. The consideration for decay of r] is 
similar. If the decay tt° —>■ 7 + X exists, one expects a flux of high energy X bosons from a 
neutrino target, since 7r° are abundantly produced in the forward direction by high energy 
( a few hundred GeV) protons in a neutrino target. If X is a long lived particle, this flux 
would penetrate the downstream shielding without significant attenuation (see section 4) and 
would be observed in a neutrino detector via the Primakoff effect, namely in the conversion 
process X — ► tt° in the external Coulomb field of a nucleus (see Figure 1). 

Because the cross section for X — >• tt° conversion is proportional to Z 2 , preferable search 
for such events is in high-Z detectors. The experimental signature of X — ► 7r° conversion is 
a single high energy 7r° decaying into two photons which results in isolated electromagnetic 
showers in the detector. The occurrence of X — * tt° conversion would appear as an excess of 
neutrino-like interactions with pure electromagnetic final states above those expected from 
SM predictions. Note that X — > 77 conversion could also be identified through charged decay 
modes of rj mesons in the final state. 

The energy spectrum of X-bosons at neutrino detector is expected to be harder than that 
of neutrinos, since there is no suppression of the high energy part due to life time of vr°(r]), 
as it is in case of charged pions those decays in flight are the main source of neutrino. Thus, 
one could expect that the highest energy bins have better sensitivity to X — > ir° conversion 
signal. 

The cross section for X — > 7r° conversion via the Primakoff mechanism is given by the 
minimal modification of the corresponding formula for the axion photoproduction fllT| |. In 
the lab system, we find 




z 



Figure 1: Feynman diagram for production of ir° by Primakoff effect. 



^{X + Z^^ + Z) = \Br{^^X 1 ) 



8r(7r° 



^f(^ do) 



where Br(ir° —>■ Xj) 



r(7r°^77) 

is the target form factor, t = — (pi — p^) 2 and A 



m^A t 2 
— ~ " A ' > , P is the momentum of the X-boson, a = Z 2 F 2 (t) 

2 

1 £) 3 . Since the minimal momentum 



square in our case t m j n = ( X 2P ^ ) 2 i is rather small, atomic form-factors must be used when 
t is small and a nuclear form factor for t > t = 7.39 • m 2 e |fTl|-[13|. It should be noted 
that in comparison with formula (18) of ref.[TT] which describes the cross section of axion 
photoproduction (7 + Z — > a + Z), a factor of ~ arises from photon and X-boson non-identity 
and the factor A in the denominator comes from the kinematics of rr — > 7 + X decay. The 
target form-factor Z 2 F 2 (t) consists of three parts. At small t < to, we use Thomas- Fermi- 
Moliere model for the atomic form-factors [[C|] : 



aH 2 



Z 2 F 2 {t) = G el {t) + G ine \t) 
G e \t) = Z 2 

G ine \t) = Z 



[l + a 2 t) 2 ' 
aft 2 



(11) 
(12) 

(13) 



where a 



111.7-Z"3 



ai 



724.2-Z~3 



(1 + aft) 2 ' 

For values of t > to, we use the elastic nuclear form- factor 
Z 2 F 2 (t) = G nud {t), (14) 
G nud (t) = (15) 



(1 + ^) 2 

where d = Q.l64:A~%GeV 2 and A is the mass number. For t- values relevant to our case the 
differential cross section % can be written in the form 

at 



d(T 1 D I 

— ~ -Brin 
dt 2 v 



SttIV^tt) Z 2 F 2 (t) 
A 7) — • ot- 



miA 



(16) 



where t ~ 2P 2 (1 — cos 8) for heavy nuclei. Using the atomic and nuclear form-factors (2-6) 
we find that the total cross section for t > t min is given by 

a (X + Z -> tt° + Z) = Br(n° -> 7 + x) 87rI V 77) . q(Gi + ^ + ^ 
Z\ aH + 1 1 1 , 



2 ft t m { n 

G 2 = |K ^ fo + l 1 ) + -5^ g-^ ], (19) 

+ 1 aito + 1 O-itmin + 1 

G 3 * Z 2 (/n(^) - |) (20) 

For heavy nuclei and for X boson masses not too close to the 7r° mass and for P < 100 
GeV, we find that a 2 £o ^ 1 5 afimin ^ 1- The inelastic atomic form-factor gives a contribution 
to the total cross section of about 1%, which we neglect. The total cross section then takes 
the simple form: 

.(* + Z + Z) - ^ X 7 )^p2) ^ • ft • ) - i] (2D 

For Pb (Z = 82, A = 207), we find a = 50.5MeV _1 , at = 75.3MeV _1 , \fd = 68.2 MeV. Note 
that the total cross section (12) depends neither on the atomic radius a, a\ nor on the value 
to (the border between the application of atomic and nuclear form-factors). 

Numerically for r(7r° — > 77) = 7.7 eV and for X boson masses Mx < 100 MeV, we find 
that the total cross section of X — > tt° conversion on lead depends rather weakly on the 
incoming X boson momentum as shown in Figure 2. 

The approximation (12) works up to Mx < 125 MeV. For Mx close to m n it is necessary 
to use formulae (8-11). The accuracy of the form-factor calculations is estimated to be better 
than 5%, [I3|; so the accuracy of our formula (12) is of the same order of magnitude. 



4 Estimate of the total cross section for X interactions 
with matter 

The total cross section a t (X + Fe — > all) for X- interactions in the Fe- shielding used in 
neutrino experiments can be estimated in the following way. We consider the particular case 
where the X boson interaction with u- and d- quarks is proportional to the electromagnetic 
interaction, namely: 

U = gxX^u^u - ^d 1>x d] (22) 
In this case, one can find that the total V-proton interaction cross section is 

a t (Xp) = l -Br{^ - 7 + X){\ - C^Lf)^a t { 1P ) (23) 

— / 77- jy 

The X boson cross section on the nuclei in the lowest order approximation is proportional 
to the atomic number in full analogy with the case of neutrino scattering (here we implicitly 
suppose that the cross section on protons is equal to the cross section on the neutrons, an 
approximation valid to within a factor of 2). 




Figure 2: Cross section for X —> ir° conversion on lead versus X boson energy calculated for 
Mx = 10 MeV and Br(ir° —>■ 7 + X) = 1. The curve is a polynomial fit to the points. 



cr t {xA) w Aa t (Xp) (24) 
Numerically for laboratory beam momenta 1 GeV< Pl < 10 3 GeV, [13] one has: 

o- t ( 7 p) = (0.12- 0.16)m6 (25) 

For Iron (A =56) and for Mx much smaller than m v we find that 

a t {xA) ss Br(7i° -> 7 + X) x 4m6 (26) 

For J Br(7r° -> 7 + X) < 10" 4 , we find 

a t (XA) < 4 x -10" 4 m6 (27) 

This estimate of the total cross section for X interaction with matter shows that the 
assumption that X bosons are penetrating particles is correct since for Br^ir — > 7 + X) < 
10~ 4 we find that the X boson mean free path in iron is > 300 km, as compared with the Fe 
and earth shielding ( total length S3 0.4 km) used for example in the CERN SPS neutrino 
beam. 



5 Rate of X — > n conversion in the NOMAD detector. 

In this section we consider the NOMAD neutrino detector at CERNjjnj as an example in 
order to estimate the rate of X — > 7r° conversion in this experiment. The neutrino beam is 
generated by 450 GeV protons delivered by the SPS to the Be neutrino target. 



The expected number of X — > 7r° events can be calculated using the following relations: 
N x __ n0 = [Br(n° -> 7 + X)] 2 • N pot ■ J f (M x )a (M X} E x )e sel dE x (28) 

where N x _ n o is the predicted number of X — > 7r° events for the given Br(7r° —>■ 7 + X), 
£ se z is the detection efficiency, fo(M x ) is the flux of X bosons per one proton on neutrino 
target (pot), N pot is the total number of pot's, a (M x , E x ) is the cross section of X — ► 7r° 
conversion on target with mass number A, calculated for Br(-n-° —>■ 7 + X)=l, and M t 
and 5 1 are fiducial mass and fiducial area of the detector target, respectively. We note that 
Br(7r° — > 7 + X) appears twice in the formula for N x _ 7T o, through the X boson flux from 
the target, and through the Primakoff mechanism. 

The flux of the X- bosons through NOMAD detector was calculated with a detailed 
GEANT [|nj simulation used to predict the neutrino flux at the NOMAD detector, and was 
found to be of the order f (M x ) — 0.1X/2.4 x 2.4 m 2 /pot for Br(n — > 7 + X)=l and with 
average X-boson energy < E x >~ 50 GeV. Note that the flux and average energy < E x > 
are weakly depended on mass M x . This results in weak dependence of limit on branching 
ratio versus M x . 

We will consider NOMAD preshower detector ]I5[ as a lead target with M t f»l ton and 
S = 2Ax 2.4m 2 . Using N pot = 10 18 pot , Br(7T° — > 7 + X) = 3x 1(T 4 , which corresponds to 
the Crystal Barrel upper limit for M x < 10MeV/c 2 , e se i = 1 and value of the cross section 
(To(Mx, Ex) = 4.0 fib taken for simplicity at E x = 5 GeV (see Figure 2) one obtains 



N x _ n o ~ J Br(7T° -> 7 + X) 2 • f ■ N pot ■ a ■ e sel ■ ■ -j- ~ 2 x 10 3 ewents/lO 18 ^ (29) 

The main contribution to the background for the X — > 7r° conversion is expected from 
the neutrino processes which have a significant electromagnetic component in the final state, 
e.g. coherent and diffractive n° production, quasi-elastic v e scattering, etc.. Using the neu- 
trino cross sections which are well known, the total number of background events N^gd was 
estimated to be ~ 10 events. Thus, if the decay 7r° — > 7 + X exists with the branching ration 
of the order of Crystal Barrel limit its effect would be easy seen at NOMAD detector. The 
limit on Br(n° — > 7 + X) could be improved by a factor of the order J N X -^/Nhkgd — 10. 
The coupling to quarks could be constrained to a x < 10~ 7 . 

The similar estimate for X — > 77 conversion rate shows that a limit of the order of the 
Crystal Barrel limit on Br(ji — ► 7 + X) could be obtained, however for a longer exposure 
time at neutrino beam, since the production of rj mesons in the forward direction in high 
energy proton collision with neutrino target is suppressed by a factor of ~ 10 2 compare to 
that of 7r°'s. 



6 Conclusion 

We have shown that a new light gauge boson X which might be produced in the decays of 
pseudoscalar mesons 7c°(r]) — > 7 + X could be effectively searched for in neutrino experi- 
ments . The X's being produced in neutrino target from these decays would penetrate the 
downstream shielding , and be observed in a neutrino detector via the Primakoff effect, the 
process of X + Z — ► tt°(i]) + Z conversion in the external Coulomb field of a nucleus. 
Our estimates show, that NOMAD experiment at CERN is able to improve current limit on 



Br(n° — > 7 + X) at least by a factor of 10 for the mass region Mx < 10 MeV during ~ 1 
month of exposure at neutrino beam. 

Note that direct searches of these rare radiative meson decays were performed in the 
previous experiments [|7| by searching for a peak in inclusive photon spectra from two-body 
decay of the tagged 7r°(77) mesons, where background due to missing of one of the two photon 
from principal decay mode vr°(^) — > 2 / y is dominated for the low Mx mass region. Our method 
is free from this disadvantage since its sensitivity is weakly depended on Mx- Note also that 
another promising way to search for a new light gauge boson from the pseudoscalar meson 
decays in neutrino experiments is to analyse the ratio R = a{vNC) /a{vCC) which might 
result in better sensitivity. However, in this case a knowledge of the final state structure in 
X interactions with a quarks is required. This subject is now under study. 
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